NT@UW-08-12 



The neutron negative central charge density: an inclusive-exclusive connection 
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Models of generalized parton distributions at zero skewness are used to relate the behavior of 
deep inelastic scattering quark distributions, evaluated at high x, to the transverse charge density 
evaluated at small distances. We obtain an interpretation of the recently obtained negative central 
charge density of the neutron. The d quarks dominate the neutron structure function for large values 
of Bjorken x, where the large momentum of the struck quark has a significant impact on determining 
the center of momentum, and thus the "center" of the nucleon in the transverse position plane. 
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Much experimental technique, effort and ingenuity has 
been used recently to measure the electromagnetic form 
factors of the nucleon [j], d, d, 0] • These quantities are 
probability amplitudes that the nucleon can absorb a 
given amount of momentum and remain in the ground 
state, and therefore should determine the nucleon charge 
and magnetization densities. 

In the non-relativistic case, the form factors are simply 
the Fourier transforms of the rest frame spatial distribu- 
tions, and the charge and magnetization mean square 
radii are derived from the slope of the form factors at 
Q 2 = 0. In the relativistic case, this interpretation is not 
correct because the wave functions of the initial and fi- 
nal nucleons have different momenta and therefore differ, 
invalidating a probability or density interpretation. This 
is addressed by working in the Breit frame, where the 
magnitude of the initial and final nucleon momenta are 
identical. However, one needs boost corrections of order 
Q 2 /m 2 , where m 2 is the mass of the constituent particles 
to which the boost is applied, to relate the rest frame 
and moving nucleon wave functions. These corrections 
are model-dependent so the use of the Breit frame 
does not provide a precise, model independent measure 
of the spatial distribution of the nucleon. 

A recent work showed that it is possible to obtain a 
model-independent nucleon charge density ((J. In the 
infinite momentum frame (IMF), the two-dimensional 
Fourier transform of the clastic form factor F\ , provides a 
model- independent transverse charge distribution, p± (b) , 
where b is the distance from the center of momentum in 
the transverse plane. The use of existing data and con- 
venient parameterizations 0, H yielded a central charge 
density of the neutron, = 0), that is negative. We 
also note that the two-dimensional Fourier transform of 
F2 can be interpreted as a magnetization density Q , and 
that this yields a difference between the magnetic and 
electric radii in the proton. 



These findings appear to contradict previous under- 
standing of the nucleon charge and magnetization dis- 
tributions based on the model-dependent extraction of 
the rest frame charge distributions. The negative core of 
the neutron transverse density also contradicts previous 
intuition that the component in which the neutron is rep- 
resented as a proton surrounded by a negatively charged 
pion cloud causes the central charge density to be posi- 
tive. This negative core is a feature even in models that 
include a pion cloud effect to reproduce the measured 
values of F". It is therefore important to understand 
the differences between this model-independent trans- 
verse charge density and the rest frame charge density 
to fully understand the new features of the transverse 
spatial distributions. 

Our goal is to obtain further information about the 
neutron charge density by using generalized parton dis- 
tributions (GPDs) which contain information about the 
longitudinal momentum fraction x as well as the trans- 
verse position b. Experimental information regarding the 
x dependence is obtained by using GPDs to reproduce 
both deep inelastic scattering and elastic scattering data. 
Thus we use this inclusive-exclusive connection to better 
understand the central neutron charge density. 

To start the analysis, we recall that form factors are 
matrix elements of the electromagnetic current operator 
J^(x u ) in units of the proton charge. The momentum 
transfer q is space-like, so that Q 2 = ~q 2 > 0. The 
normalization is such that -Fi(O) is the nucleon charge, 
and ^2(0) is the proton anomalous magnetic moment. 
The Sachs electric and magnetic form factors are given 
by G E = F 1 - (Q 2 /4M 2 )F 2 and G M = F 1 +F 2 . 

The widely studied GPDs [H 03 are of high current 
interest because they can be related to the total angular 
momentum carried by quarks in the nucleon. We con- 
sider the specific case in which the longitudinal momen- 
tum transfer £ is zero, and the initial and final nucleon 
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hclicities are identical (A' = A). Then, in the light-cone 
gauge, A + = 0, the matrix element defining the GPD H q 
for a quark of flavor q and zero skewness is 

= |^ r (p + , P ',A|6 g ( a; ,0)b + ,p,A)e^ + ^,(l) 

where 

d 9 (x,b) = | ^ (-T' b ) 9+ (t ,1> ) e " P+ "-( 2 ) 

We abbreviate H q (x, £=0,t) = H q (x,t) and —t = — (p' — 
P) 2 = (p'~p) 2 = Q 2 ■ The simple form of t results from its 
invariance under transverse boosts [l2j : Lorentz transfor- 
mations, defined by a transverse vector v that transform 
a four- vector k according to k + — > fc + , k — > k /c + v 
and fc _ such that k 2 is unchanged. These quantities are 
part of a kinematic subgroup of the Poincare group that 
obey the same commutation relations as those among 
the generators of the Galilean transformations for non- 
relativistic quantum mechanics in the transverse plane. 

GPDs allow for a unified description of a number 
of hadronic properties [Toj . The most relevant for us 
arc that for t=0 they reduce to conventional PDFs, 
H q (x,0) = q(x), and that the integration of the charge- 
weighted H q over x yields the nucleon electromagnetic 
form factor: 

F 1 (t) = J2e q I dxH q (x,t). (3) 
g J 

The s pat ial structure of a nucleon can be exam- 
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using nucleonic states that are trans- 
versely localized. The state with transverse center of 
mass R set to 0, |p + ,R = 0, A) is formed by taking a 
linear superposition of states of transverse momentum. 
Doing this requires the use of a frame with infinitely large 
P + - 

The impact parameter-dependent PDF [16[ is the 
matrix clement of the operator O q in the state 
|p+,R = 0,A): 

p\{b,x) = (p + ,R = 0,A|6 f; (x,b)|p+,R = 0,A). (4) 

We use the notation p q ± (h, x) instead of the originally de- 
fined [l6j q(x, b) because the quantity truly is a density, 
giving the probability that the quark has a longitudi- 
nal momentum fraction x and is at a transverse position 
b. The quantity p q ± (h,x) is the two-dimensional Fourier 
transform of the GPD H q : 

Pi(b,x) = J ^e- i ^ h H q (x,t = - q 2 ), (5) 

with H q appearing because the initial and final hclicities 
are each A. A complete determination of H q (x,t) (with 
t < 0) would determine p±(x,h). 

One finds can extract the form factor Fi [3] by inte- 
grating p^(b, x) over all values of x, multiplying by the 



quark charge e g , and summing over quark flavors q. The 
resulting IMF charge density in transverse space is 

pf(b) = £e, f dxpl(b,x) 

This quantity gives the charge density at a transverse 
position b irrespective of the longitudinal momentum 
fraction. The primary difference between the present 
charge density p±{b) and the older interpretation that 
the charge density as the three-dimensional Fourier trans- 
form of Ge is that the present approach provides a 
model- independent, two-dimensional charge distribution 
in the plane transverse to the motion of the nucleon in the 
infinite momentum frame. The boost corrections here are 
simply kinematic and are incorporated in the formalism. 
In the older interpretation, the model-dependent boost 
corrections can not be avoided. 

Our aim is to investigate px(h, x) to understand the 
origin of the neutron's negative central charge den- 
sity. The quantities are not measured directly, but have 
been obtained from models that incorporate fits to par- 
ton distributions and electromagnetic nucleon form fac- 
tors [l?], [H, O H3] • This method exploits form factor 
sum rules at zero skewness, obtained neglecting the ef- 
fects of strangeness, to obtain information regarding the 
valence quark GPDs, H« = - H*. This yields the net 
contribution to the form factors from quarks and anti- 
quarks, although it does not correspond to the valence 
distribution within a model for which sea distributions 
for quarks and antiquarks have different i or i depen- 
dences. To proceed further one must model the GPDs, 
and the results can be expected to depend on the chosen 
forms. Diehl et al. [ljj use 

H%{x,t)=q v {x)e X p[f q (x)t], (7) 

where 

f q (x) = [a'log[l/x] + B q ](l - x) 3 + A q x(l - x)\ (8) 

is the form that gives the best fit to the data. The pa- 
rameter a' represents the slope of the Re gge trajectory 
{a' = 0.9 GeV 2 ), and the CTEQ6 PDFs [2l| are taken 
as input. Here we use the best fit p arameters, taken 
from the second line of Table 8 of [17| ■ These are A u = 
1.26 GcV~ 2 : B u = 0.59 GcV~ 2 , A; = 3.82 GcV~ 2 ,B d = 
0.32 GeV -2 . We note that the labels u and d here re- 
fer to the u and d quarks in the proton. These corre- 
spond to d and u quarks in the neutron, if charge sym- 
metry [H, [H, [H, 111] is upheld. It is well-known that 
for the proton, 2d v /u v falls rapidly for large values of x, 
which means that u quarks dominate the parton distribu- 
tion for large values of x. This means that in the neutron, 
the d quarks dominate the parton distribution for large 
values of x. The distributions of [l8[ have A„ = B q = 
and f q (x) = [a' q log[l/x]](l — x). Those of [19| have a 
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more complicated form and also include the constraint 
that the nucleoli consists of three quarks at an initial 
scale of Ql = 0.094 GeV 2 . 
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FIG. 1: (Color Online) The proton transverse charge density, 
p^(6, Arr), for quarks in different Ax regions: a;<0.15 (solid), 
0.15<x<0.3 (long-dash), 0.3<x<0.5 (short-dash), and x>0.5 
(dotted). The curves are calculated from the GPD of Ref. [l7l | 
and have been normalized to unity at b = 0. 



Now consider the charge distribution of the neutron. 
We expect that the d quarks dominate at large x and 
therefore become important at small values of b. Because 
the distribution of quarks at large x will be highly local- 
ized near b = 0, a negative peak can be formed if the large 
x distribution is sufficiently dominated by down quarks, 
thus yielding a significant contribution of negative charge 
at large enough x. At very low x values, the valence dis- 
tribution for up quarks in the neutron is roughly half 
that of the down quarks, d™(x)/u™(x) « 2, and the net 
charge coming from u and d quarks will approximately 
cancel, although the distribution as a function of b need 
not be zero everywhere. Above x = 0.5, d"(x) is at least 
three times the size of u™(x), and increases with x. So 
for x > 0.5 the net impact to the charge distribution will 
be negative, and will be peaked at smaller values of b. 
We show this explicitly in Fig. [2l where we separate the 
contributions to the neutron charge density from u and 
d quarks based on the GPD fit of Ref. 17| . The distri- 
butions of and [l9| yield somewhat different results, 
but they exhibit the same qualitative behavior. 



Our goal here is to examine the connection between 
regions of x and regions of b. To do this we define 



P 9 ±(b,Ax) 



dx e q p\{b,x), 



(9) 



with pj" being obtained from appropriate sums of p q j_. 
This represents the contribution to the charge density 
from quarks in the x region defined by Ax. 

An important feature of the present approach is that 
these charge distributions are taken with respect to the 
center of momentum in the transverse plane. Thus, 
the transverse position b is taken with respect to the 
momentum-weighted average position of all partons, in- 
cluding the struck quark. At low x, the struck quark 
has little impact on the center of momentum, and this 
corresponds to intuitive picture of spatial distribution. 
At large x, the struck quark plays a significant role in 
defining the CM, and so distribution becomes localized 
at small values of b. This can be seen in Fig. [T] where 
for x ~ 0.1, the half-maximum width is 0.5 fm, while for 
x w 0.8, it is 0.12 fm. The curves have been scaled to 
yield unity at b = 0, to emphasize the variation in width. 
The four Ax regions yield 58%, 25%, 14%, and 3% of the 
total charge, with the largest contributions coming from 
the bins with the smallest values of x. Thus the large 
x quarks, dominantly u quarks in the proton, play an 
increasingly prominent role in the charge distribution at 
small values of b. The figure obtained using the Guidal et 
al, parameterization for the GPDs is barely distinguish- 
able from Fig. □ The GPDs of [3 also have a strong 
tendency to be constrained to smaller and smaller values 
of b as the value of x increases. We evaluate the GPDs 
of all three models using the starting scale Qq of each 
model. 
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FIG. 2: (Color Online) The u and d quark contributions to the 
neutron transverse charge density, p^(b,Ax) and pi (6, Ax). 
The curves correspond to the same Ax regions as in Fig. [T] 
The largest contributions come from small x, where u and d 
quarks contribute roughly equal amounts of charge. As one 
goes to larger x values, the charge is shifted to smaller values 
of b, while at the same time the contribution from the up 
quarks drops rapidly with respect to the down quarks, due to 
the rapid falloff of the neutron u to d quark ratio at large x. 



The next step is to examine the total charge distri- 
bution of the neutron. Fig. [3] separates the contributions 
from low and high x regions. For x < 0.23 the charge dis- 
tribution is positive for b < 1.5 fm and slightly negative 
distribution at larger radii. For x > 0.23, the contribu- 
tion is largely negative, and highly localized below 0.5 fm. 
The negative region at the center of the neutron trans- 
verse charge distribution arises a natural consequence of 
the model-independent definition of the charge density. 
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The low momentum partons have a larger spatial extent 
and reproduce the intuitive result of the pion cloud pic- 
ture: a positive core with a small negative tail at large 
distances. 
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FIG. 3: (Color online) Transverse charge density for the neu- 
tron. The dotted line is the contribution from x < 0.23, 
dashed is that for x > 0.23, and the solid is the total. 



A more intuitive picture of the charge distribution can 
be obtained by looking at the distribution of charge rel- 
ative to the spectator partons, so that the struck quark 
does not influence the definition of the center of mass. 
This can be approximated by looking at the position of 
the struck quark relative to the spectators. We work in 
the transverse plane, with the origin set to the center 
of momentum, giving a^bi = 0. For a struck quark 
at (xi,bi) = (x, b), we can determine the momentum- 
weighted spectator position, b spC c, and the relative dis- 
tance from the struck quark to the spectator quarks: 

xibi + ^ x i^i =xh+(l- z)b S pcc = 0, (10) 

i>l 

B re ] = b b spcc = — = B re i. (11) 

(1-x) 

We exhibit the dependence on B rG i by defining a function 
^ poc (B. rei ,x) =p x (B rel {l-x),x) (12) 

which gives the probability that a struck quark of lon- 
gitudinal momentum fraction x is a distance B re i away 
from the spectator center of momentum. Figure [4] shows 
this rescaled version of p_i_(6), with the contribution at 
each x value normalized to unity at b = 0. The quantity 
p£>P ec (B ro i, x ^ can no ^ determined in a model indepen- 
dent manner, but may be a better approximation to our 
intuitive picture of the charge distribution, as it removes 
the influence of the struck quark on defining the center of 
the nucleoli. While the charge distribution coming from 
very low x quarks has a greater spatial extent, the de- 
creasing width of the p^ (b) distribution for large x quarks 
is essentially completely removed when looking at -B re i- 



Before concluding, it is worthwhile to comment on the 
relation between the present work and the difference be- 
tween the electric and magnetic radii of the proton [§]. 
In the model- independent, IMF approach presented here, 
the electric and magnetic transverse radii have a clear 
connection to F\ and F2 and a Foldy [2(| term causes a 
difference between the transverse radii. The Foldy term 
is responsible for most of the charge radius defined by 
Ge- Understanding the neutron's negative central den- 
sity is more subtle and requires knowledge of p(x, b). 




FIG. 4: (Color online) The u and d quark contributions to 
p^ pec '™(S re i,x) see Eq. dUJ). vs B ro i for x = 0.1 (solid), 
0.3 (long-dash), 0.5 (short-dash), and 0.7 (dotted). The 
curves are scaled to unity at B IC \ — 0. 



We summarize our findings with the statement that, 
using the model GPDs of Refs. [H GJ, El, the domi- 
nance of the neutron's d quarks at high values of x leads 
to a negative contribution to the charge density which, 
due to the definition of b, becomes localized near the 
center of mass of the neutron. This localization does not 
appear when examined as a function of the position of 
the struck quark relative to the spectators. 



Acknowledgments 

This work was supported by the U. S. Department 
of Energy, Office of Nuclear Physics, under contracts 
FG02-97ER41014 and DE-AC02-06CH11357. We thank 
D. Geesaman, P. Kroll and B. Wojtsekhowski for useful 
discussions. We thank the ECT* for hosting a work- 
shop where many of the calculations we present were 
performed. 



[1] H. Gao, Int. J. Mod. Phys. E12, 1 (2003 [Erratum-ibid., [14 

567, 2003]). [15 
[2] C. E. Hyde-Wright and K. de Jager, Ann. Rev. Nucl. [16 

Part. Sci. 54, 217 (2004). [17 
[3] C. F. Perdrisat, V. Punjabi, and M. Vanderhaeghen, 

Prog. Part. Nucl. Phys. 59, 694 (2007). [18 
[4] J. Arrington, C. D. Roberts, and J. M. Zanotti, J. Phys. 

G34, S23 (2007). [19 
[5] J. J. Kelly, Phys. Rev. C 66, 065203 (2002). [20 
[6] G. A. Miller, Phys. Rev. Lett. 99, 112001 (2007). 
[7] J. J. Kelly, Phys. Rev. C 70, 068202 (2004). [21 
[8] R. Bradford, A. Bodek, H. Budd, and J. Arrington, Nucl. [22 

Phys. Proc. Suppl. 159, 127 (2006). 
[9] G. A. Miller, E. Piasetzky, and G. Ron (2007), 0711.0972. [23 
[10] X.-D. Ji, Phys. Rev. D55, 7114 (1997), hep-ph/9609381. [24 
[11] A. V. Radyushkin, Phys. Rev. D56, 5524 (1997), hep- 

ph/9704207. [25 
[12] J. B. Kogut and D. E. Soper, Phys. Rev. D 1, 2901 

(1970). [26 
[13] M. Diehl, Eur. Phys. J. C25, 223 (2002). 



D. E. Soper, Phys. Rev. D 15, 1141 (1977). 

M. Burkardt, Int. J. Mod. Phys. A18, 173 (2003). 

M. Burkardt, Phys. Rev. D 62, 071503(R) (2000). 

M. Diehl, R. Feldmann, Th. an Jakob, and P. Kroll, Eur. 

Phys. J. C39, 1 (2005). 

M. Guidal, M. V. Polyakov, A. V. Radyushkin, and 

M. Vanderhaeghen, Phys. Rev. D 72, 054013 (2005). 

S. Ahmad et al., Phys. Rev. D 75, 094003 (2007). 

B. C. Tiburzi, W. Detmold, and G. A. Miller, Phys. Rev. 

D70, 093008 (2004), hep-ph/0408365. 

J. Pumplin et al., JHEP 07, 012 (2002). 

G. A. Miller, B. M. K. Nefkens, and I. Slaus, Physics 

Reports 194, 1 (1990). 

G. A. Miller, Phys. Rev. C 57, 1492 (1998). 

J. T. Londergan and A. W. Thomas, Prog. Part. Nucl. 

Phys. 41, 49 (1998), hep-ph/9806510. 

G. A. Miller, A. K. Opper, and E. J. Stephenson, Ann. 

Rev. Nucl. Part. Sci. 56, 253 (2006). 

L. L. Foldy, Phys. Rev. 83, 688 (1951). 



